AO-AIU  79%  NAVAL  SURFACE  WEAPONS  CENTER  SILVER  SPR1N6  MD  F/6  20/4  \ 

COMPARISON  OF  NUMERICAL  RESULTS  ANO  MEASURED  DATA  FOR  SMOOTH  AN— €TC<U) 

MAY  61  T  HSIEH 


UNCLASSIFIED 


NL 


ac 


COMPARISON  OF  NUMERICAL  RESULTS  AND  MEASURED  DATA  FOR 
SMOOTH  AND  INDENTED  NOSETIPS 


Tsuying  Hsieh 

Naval  Surface  Weapons  Center,  White  Oak,  Maryland 


ABSTRACT 


/f 


Numerical  calculations,  using  an  unsteady  implicit  numerical 
algorithm  which  solves  either  the  inviscid  or  the  thin-layer 
Navier-Stokes  equations,  were  performed  for  smooth  and  severely 
Indented  nosetips  at  hypersonic  speed  and  zero  angle  of  attack. 

The  computed  results  of  inviscid  and  laminar  flowfield  are 
compared  to  wind  tunnel  measured  data  for  surface  pressure,  shock 
location,  heat  transfer  and  density  distribution  in  the  shock  layer. 
Good  agreements  between  the  calculated  and  measured  flowfield  are 
obtained  for  smooth  nosetlp  without  flow  separation.  Difficulties 
in  the  simulation  of  severely  indented  nosetips  with  large  separation 
bubble  or  sharp  corner  are  discussed.  ^ 

NOMENCLATURE  ... 
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cp  -  Pressure  coefficient,  p-p^fh  P*,  U*,  2 

J,K  =  Index  of  grid  in  the  streamwise  and  .. 
body-to-shock  direction  respectively 

L  *  Reference  length,  2.56  cm  (1  in.) 

M  ■  Mach  number 

n  *  normal  direction 
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U  ■  Velocity 


X,Y  *  Axial  and  radial  coordinates  respect¬ 
ively 

p  *  Density 

eE  ■  Coefficient  for  explicit  dissipation 
eI  ■  Coefficient  for  implicit  dissipation 
V  ■  Kinematic  viscosity  of  fluid 
Subscript 
<*>  ■  free  stream 
l  «  laminar  flow 
o  *  stagnation  point 
t  *  turbulent  flow 
w  *  wall 


INTRODUCTION 


Because  of  ablation,  the  nosetip  of  a  spherical  body  undergoes  continuous 
change  during  re-entry.  The  shape  of  the  nosetip  has  a  great  influence  on 
the  flowfield  over  the  entire  body,  i.e.,  the  nose  region  and  thus  the 
afterbody.  In  order  to  understand  the  flowfield  about  indented  nosetip  shapes 
that  are  likely  to  occur  during  the  re-entry  process,  considerable  effort  has 
been  expended  both  experimentally  and  theoretically. (also  see  references 
listed  in  Ref.  1  and  2). 


Among  the  experimental  work.  Refs.  3  and  4  provide  a  comprehensive  set  of 
measurements  for  a  systematic  change  of  nosetip  shapes.  This  Includes 
surface  pressures,  heating  rates,  flowfield  pictures  using  electro-optical 
techniques  and  velocity  mapping  using  Laser  Doppler  Velocimeter.  These 
experiments  are  intended  to  provide  basic  information  about  the  important  flow 
features  to  be  expected  in  the  flowfield  and  also  serve  as  a  useful  data 
base  for  the  development  of  a  numerical  code  to  predict  the  flowfield. 

Among  the  many  numerical  schemes  intended  for  Indented  nosetip  calculation, 
a  promising  and  versatile  one  seems  to  be  due  to  Kutler  et  al . *  As  described 
in  Ref.  1,  Kutler  et  al  solved  the  unsteady  Navler-Stokes  equations  with  the 
"thin-layer"  approximation  for  nosetip  of  arbitrary  shapes  at  zero  lndlcence 
using  the  implicit  factored  numerical  algorithm  of  Warming  and  Beam. 5  The 
steady  solution  is  obtained  asymptotically  in  time  and  both  viscous  and 
invlscld  flowfields  can  be  computed  using  the  same  computer  program,  refered 
toasK-C-L  code  in  this  paper.  As  described  in  Ref.  2,  when  the  inviscid  - 
portion  of  the  K-C-L  code  was  applied  to  the  nosetip  shapes  reported  in 
Ref.  3  and  4  surmountable  difficulties  were  encountered  during  the  course  of 
calculation  because  of  the  presence  of  small  radius  expansion  corners  and  a 


concave  compression  turn  In  these  nosetips.  It  was  later  found  that  a  special 
calculation  procedure  is  required  in  order  to  obtain  reasonable  solutions. 

For  the  details  of  the  numerical  procedure  used  in  the  K-C-L  code  and  the 
special  calculation  procedure  for  the  indented  nosetip  computation,  the 
readers  should  consult  to  Ref.  1  and  2  and  will  not  be  repeated  in  this  paper. 
The  same  calculation  procedure  has  since  been  applied  to  compute  viscous  flow 
over  smooth  and  Indented  nosetips  and  several  typical  results  will  be  described 
and  compared  to  measured  data. 


II.  RESULTS  AND  DISCUSSION 


A.  Smooth  Nosetips 

The  comparison  of  inviscid  results  for  surface  pressure  and  density 
distribution  over  a  sphere  is  reproduced  here  from  Ref.  2  as  shown  in  Fig.  1. 

The  grid  used  in  the  computation  is  J  *  32  (along  body  surface)  x  K  *  12 
(across  the  shock  layer)  and  the  total  number  of  time  steps  used  is  600 
(with  a  courant  number  of  2.5).  This  requires  about  2  minutes  CPU  time  in 
CDC-7600  computer.  The  surface  pressure  is  compared  to  the  work  of  Inouye 
and  Lomax<>  and  the  measured  data  of  Baer?  as  shown  in  Fig.  la.  The  density 
distribution  in  the  shock  layer  is  compared  to  the  measured  data  of 
Sedney  and  Kahl®  as  shown  in  Fig.  lb.  The  agreement  are  seen  to  be  very  good. 

Viscous  flow  calculations  based  on  the  thin-layer  approximation  of  Navier-Stokes 
equations,  were  performed  for  a  hemisphere-cone  under  the  following  flow 
condition:  Moo  =  5.92,  Re*,  *  5.5  x  106/ft  and  Tw/Ta>  *  4.78.  A  grid  of  28x32 
was  used.  For  laminar  calculation,  400  time  steps  with  a  Courant  number  (CN) 
of  75  can  be  considered  as  the  final  solution  (the  nondimensional  shock  speed 
is  in  the  order  of  10“^  and  the  CPU  time  Is  about  5  min).  As  shown  in  Fig.  2 
the  calculated  results  for  heat  transfer  in  term  of  Stanton  number  over  the 
surface  is  compared  to  the  measured  data  reported  in  Ref.  7.  Also  plotted 
in  Fig.  2  is  the  boundary  layer  calculation  using  Cebeci-Smith’ s  boundary 
layer  code  as  given  in  Ref.  7.  It  is  seen  that  the  agreement  is  satisfactory. 

It  should  be  noted  that  for  hemisphere-cone,  the  flow  is  fully  attached.  The 
surface  pressure  agrees  well  among  the  inviscid  and  laminar  calculation  as 
well  as  the  measured  data.  Therefore,  the  good  agreement  in  heat  transfer 
between  the  boundary  layer  and  the  thin  layer  N-S  calculation  is  to  be  expected. 

In  the  numerical  procedure,  there  are  two  kinds  of  dissipation  used  for 
stability  purpose,  l.e.,  the  explicit  and  Implicit  dissipation  terms  with 
coefficients  of  CE  and  EI  respectively.  For  inviscid  and  laminar  calculation, 
solution  can  be  obtained  without  the  implicit  dissipation.  For  turbulent  calculation, 
it  was  found  that  the  implicit  dissipation  must  be  added  in  order  to  have  a 
converged  solution.  Since  the  Implicit  dissipation  terms  have  no  effect  on  the 
solution  as  the  steady  state  is  approached,  it  is  therefore  used  in  both  laminar  and 
turbulent  calculation.  The  explicit  dissipation  terms,  on  the  other  hand, 
will  always  present.  In  viscous  flow  calculation,  it  is  desirable  that  the 
explicit  dissipation  coefficient  be  kept  minimum  so  that  the  real  viscous 
effects  can  be  simulated.  From  the  experience  for  indented  nosetip  calculation, 
the  value  of  eE  cannot  be  too  small  in  order  to  have  a  smooth  bow  shock  (no 
wiggling)  where  the  flow  can  essentially  be  assumed  to  be  inviscid.  Therefore, 
the  eE  value  in  the  n  direction  is  linearly  reduced  to  one  tenth  of  the  input 
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value  toward  the  wall  in  order  to  minimize  its  effects  there.  This  devise 
also  consistent  with  the  thin-layer  approximation.  When  eE  is  linearly 
reduced,  it  is  denoted  by  £E  '  to  distinguish  from  the  uniform  one.  As  shown 
in  Fig.  2,  the  effects  of  eE  '  are  seen  to  be  insignificant  in  the  heat  transfer 
results  for  hemisphere-cone  with  eE  "  varies  from  0.1  to  0.02  and  GI  *  3  eE '. 
This  is  not  so  for  indented  nosetip  calculation. 

In  the  simulation  of  turbulent  flow,  the  turbulence  model  of  Baldwin  and 
Lomax ^  was  used.  The  final  laminar  solution  with  eE  *  0.1  was  used  as  the 
initial  flowfield  for  turbulent  calculation.  A  turbulent  solution  was 
produced  in  200  time  steps  with  the  nondimen sional  shock  speed  in  the  order 
of  10~3.  As  shown  in  Fig.  2,  the  Stanton  number  increase  significantly  for 
turbulent  flow  as  compared  to  the  laminar  flow  one  (no  measured  data  for 
comparison).  It  should  be  pointed  out  that  the  surface  pressure  obtained  from 
the  laminar  and  turbulent  solution  agree  to  two  digits. 


B.  Indented  Nosetips 

The  special  calculation  procedure  for  inviscid  flow  over  indented  nosetips 
as  described  in  Ref.  2  was  applied  to  calculate  viscous  flow  over  indented 
nosetips.  Model  1  and  4  as  described  in  Ref.  2  were  chosen  for  this  investigation. 
The  simulation  of  viscous  turbulent  flow  over  indented  nosetips  has  not  been 
very  successful  as  described  in  the  following  paragraphs. 

Laminar  flow  over  Model  4  was  first  calculated.  The  calculation  started 
with  a  grid  of  24  x  32  (CN  ■  150,  eE  *  0.4,  eI  *  0)  for  400  time  steps  to 
obtain  a  laminar  solution  over  a  sphere  at  M«  *  5.0,  Re  ■  8  x  106/FT  and 
Tw/T„  ■  5.4.  The  sphere  was  then  deforned  to  the  shape  of  Model  4  in  1800 
time  steps.  The  grid  was  then  increased  to:  (!)  58  x  32  and  (B)  56  x  48  in 

another  1600  time  steps  each  with  the  final  values  of  CN  *  50,  £E  ■  0.1  for 
(A)  and  eE  '  •  0.3  for  (B).  The  calculated  surface  pressure  and  shock 
locations  from  these  two  solutions  are  close  (results  not  shown).  This  provides 
a  self  verification  of  the  results.  Since  grid  (B)  contains  more  points  in 
the  n  direction,  its  solution  resolves  the  flow  field  better  and  is  used  for 
comparison  as  shown  in  Fig.  3  and  4. 

Unlike  the  hemisphere-cone,  the  turbulent  calculations  for  Model  4 
encountered  serious  difficulties.  Large  amplitude  oscillation  of  pressure 
in  the  flowfield  quickly  interrupted  the  computation.  The  value  of  CN  was 
gradually  reduced  and  the  value  of  eI  was  increased  (eE  '  *  0.3  was  maintained). 

At  CN  -  2  and  CI  -  6,  it  was  possible  to  run  for  200  time  steps  with  the 
non-dimensional  shock  speed  coverging  to  a  value  of  0.04.  The  shock  speed  is 
then  started  to  Increase  slowly  but  steadily.  A  further  Increase  of  GI  up  to 
12  would  not  help  to  obtain  a  converged  solution.  Thus,  the  solution  before 
the  shock  speed  started  to  increase  is  shown  in  Fig.  3  and  4  for  comparison. 

It  is  understood  that  the  turbulent  solution  shown  in  Fig.  3  and  4  is  not 
the  converged  solution. 

Fig.  3  shows  the  comparison  of  shock  location  between  calculations  and 
experiments.  The  inviscid  shock  layer  is  thinner  around  the  indented  region 
as  expected.  The  laminar  and  turbulent  solutions  for  shock  location  are 
almost  coincide  and  fall  in  between  the  inviscid  solution  and  the  measured 
data.  The  primary  separation  bubble  indicated  by  the  laminar  solution  is 


smaller  than  observed  experimentally.  The  laminar  separation  point  of  the 
primary  separation  bubble  is  at  the  downstream  end  of  the  expansion  corner, 
but  the  flow  picture  (Fig.  8a  of  Ref.  2)  shows  that  the  flow  separates 
immediately  at  the  beginning  of  the  expansion  corner.  In  general,  the 
effects  of  turbulence  is  to  move  the  separation  point  toward  downstream  and 
the  separation  bubble  will  be  smaller.  This  fact  suggests  that  the  discrepancy 
is  not  because  of  turbulence  effects  but  from  some  other  sources  which  are  not 
captured  in  the  numerical  solution.  Also  the  laminar  solution  indicates  that 
there  is  a  secondary  separation  bubble  within  the  primary  separation  bubble 
as  sketched  in  Fig.  3.  The  secondary  separation  bubble  was  not  reported 
in  Ref.  4. 

In  Fig.  4,  the  surface  pressure  distribution  obtained  from  the  inviscid, 
laminar  and  turbulent  solution  are  compared  to  the  measured  data.  It  is 
noted  that  the  viscous  solutions  compare  better  with  the  measured  data  than 
the  inviscid  curve.  The  region  around  the  expansion  corner  S/L  0.26  -  0.7 
where  the  Inviscid  and  viscous  solutions  are  seen  to  agree  well  (i.e.,  no 
flow  separation)  but  are  lower  than  the  measured  data.  The  dip  in  the 
pressure  curve  in  the  region  S/L  'u  2.5  (where  the  secondary  separation 
bubble  starts)  is  not  shown  by  the  measured  data.  A  calculation  was  made  with 
the  complete  Navier-Stokes  equation  for  laminar  flow  with  the  same  grid 
distribution  as  used  in  the  thin-layer  approximation  calculation  and  the 
results  for  surface  pressure  agree  up  to  two  digits. 


For  Model  1,  there  is  a  sharp  expansion  corner  with  a  radius  of  0.062  inch. 

Three  grid  points  were  used  to  cover  the  corner  as  was  done  for  the  inviscid 
calculation  reported  in  Ref.  2  and  a  total  grid  points  of  33x32  were  used  for  the 
viscous  calculation.  Only  laminar  solution  can  be  obtained.  As  shown  in  Fig  5 
and  6,  the  inviscid  solution  agrees  reasonably  well  with  the  measured  data  for 
both  the  shock  location  and  the  surface  pressure  because  the  separation  bubble 
is  small  as  indicated  by  the  flow  picture  (Fig.  5a  of  Ref.  2).  In  contrast  to 
the  result  of  Model  4,  however,  the  laminar  solution  indicates  that  the  flow 
separates  immediately  after  the  comer  and  form  a  large  primary  separation  bubble 
as  shown  in  Fig.  5.  Within  the  primary  separation  bubble,  there  is  also  a 
secondary  separation  bubble  around  the  location  of  compression  turn.  As  a  result 
of  the  primary  separation  bubble,  the  laminar  shock  layer  becomes  thicker  near  the 
separation  bubble  and  thinner  afterwards  as  compared  to  the  measured  data.  Also  the 
surface  pressure  obtained  from  the  laminar  solution  looks  entirely  wrong. 

It  was  not  possible  to  obtain  a  turbulent  solution  for  Model  1,  not  even 
one  like  that  of  Model  4.  The  obvious  reason  is  that  the  laminar  solution  is 
too  far  off  from  the  measured  data,  which  is  assumed  to  be  close  to  the  turbulent 
solution,  therefore  the  starting  flowfield  is  too  poor  to  carry  through  the 
calculation. 
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From  the  above  comparison  of  results,  it  is  concluded  that  the  numerical 
code  works  well  for  smooth  nosetips.  For  those  indented  nosetips  investigated 
in  this  paper,  the  code  fail  to  give  a  satisfactory  solution.  Possible  reasons 
for  the  failure  of  the  code  to  obtain  a  solution  and  for  the  discrepancy 
between  the  viscous  calculation  and  the  experimental  data  are  suggested  as 
follows:  (1)  The  turbulence  model  used  is  Inadequate.  (2)  The  first- 
order-accurate  numerical  procedure  is  not  sufficient  to  simulate  the  complicated 
flowfleld.  (3)  The  grid  points  and  their  distribution  are  insufficient  to 
resolve  the  viscous  effects.  (4)  The  sharp  corner  as  described  in  Model  1 
requires  special  treatment  to  avoid  flow  separation  introduced  from  sources 
other  than  the  viscous  effects.  (5)  Because  of  the  prevailing  of  a  large 
separation  bubble,  the  thin-layer  approximation  may  lead  to  wrong  solution, 

i.e.,  for  a  significant  portion  of  the  flowfleld,  the  viscous  terms  in 
both  the  £  and  n  directions  should  be  kept. 


REFERENCES 

1.  Kutler,  P.,  Chakravarthy ,  S.  R.  and  C.  P.  Lombard  "Supersonic  Flow  Over 
Ablated  Nosetips  Using  an  Unsteady  Implicit  Numerical  Procedure" 

AIAA  Paper  78-213,  Jan  16-18,  1978. 

2.  Hsieh,  T. ,  "Numerical  Investigation  of  Flowfleld  About  a  Series  of 
Indented  Nosetips,"  AIAA  Paper  No.  81-0077,  Jan.  1981. 

3.  Ragsdale,  W.  C.  and  Morrison,  A.  M.,  "IAP  202  Heat  Transfer  and  Pressure 
Tests  in  the  NSWC/WOL  Hypersonic  Tunnel,"  NSWC/WOL  MP  78-18,  October  27,  1978. 

4.  Tanta,  W.  J.,  "Indented  Nose  Flowfleld  Tests,"  WTR  1329,  Naval  Surface 
Weapons  Center,  White  Oak,  MD.,  July  1980. 

5.  Warming,  K.  F.  and  Beam,  R. ,  "On  the  Construction  and  Application  of 
Implicit  Factored  Schemes  for  Conservation  Laws,"  SIAM,  ATNS  Proceedings, 

Vol.  11,  1978,  pp.  85-129. 

6.  Inouye,  M.  and  Lomax,  H.,  "Comparison  of  Experimental  and  Numerical 
Results  for  the  Flow  of  a  Perfect  Gas  About  Blunt-Nosed  Bodies,"  NASA  TN 
D-1426,  Sept  1962. 

7.  Baer,  A.  L.,  "Pressure  Distribution  on  a  Hemisphere  Cylinder  at  Supersonic 
and  Hypersonic  Mach  Numbers,"  AEDC  TN-61-96,  Arnold  Engineering 
Development  Center,  1961. 

8.  Sedney,  R.  and  Kahl,  G.  D.,  "Interferometric  Study  of  the  Blunt  Body 
Problem,"  Ballistic  Research  Laboratory  Report  No.  1100,  1960. 

9.  Baldwin,  B.  S.  and  Lomax,  H.  "Thin  Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flow"  AIAA  Paper  78-257,  1978. 


11-302 


Comparison  of  Surface  Pressure  and  Density  Field 
Between  Calculation  and  Experiment 


5.92,  R^oo  =  5.5  x  106/FT,  TW/T  =  4.78 


FIGURE  2  COMPARISON  OF  HEAT  TRANSFER  MEASUREMENTS  ON  A  HEMISPHERE  CONE 


Figure  3  Comparison  of  Shock  Location  Between  Calculation 
and  Experiment  for  Model  4 


Re  =  8  x  106/FT,  Tw/T 


Figure  A  Comparison  of  Surface  Pressure  Between  Calculation 
and  Experiment  for  Model  A 


.  Tw/T  =  4.78 


